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Abstract 



We address the Gross-Pitaevskii (GP) equation with a periodic linear potential 
and a periodic sign- varying nonlinearity coefficient. Contrary to the claims in the 
previous works of Abdullaev et al. [PRE 77, 016604 (2008)] and Smerzi & Trom- 
bettoni [PRA 68, 023613 (2003)], we show that the intersite cubic nonlinear terms 
in the discrete nonlinear Schrodinger (DNLS) equation appear beyond the applica- 
. bility of assumptions of the tight-binding approximation. Instead of these terms, 

' for an even linear potential and an odd nonlinearity coefficient, the DNLS equation 

\^ . and other reduced equations for the semi-infinite gap have the quintic nonlinear 

f^) ' term, which correctly describes bifurcation of gap solitons. 

on" 

Keywords: Gross-Pitaevskii equation, discrete nonlinear Schrodinger equation, 
^ap solitons, bifurcations, semi-classical limit, Wannier functions. 

■ 1 Introduction 

The generalized DNLS (discrete nonlinear Schrodinger) equation with intersite cubic 
nonlinear terms, 

ic n = a(c n+1 + c n _i) + P\c n \ 2 c n 

+7( 2 l c n| 2 (c„+l + C n -i) + Cn(c n+ i + C„_i) + \c n+ i\ 2 C n+ i + |c„_i | 2 C„_i) 
> 



+S((c 2 n+1 + cl_ x )c n + 2(|c n+ i| 2 + |c n _i| 2 )c n ), (1) 



where (a, f3, 7, 5) are constant parameters and the dot denotes differentiation in time, 
was derived independently in various contents. Smerzi & Trombettoni [22] suggested 
that this equation models Bose-Einstein condensates in a lattice, when Wannier func- 
tions associated with a periodic potential are replaced by the nonlinear bound states. 
Independently, this equation was derived heuristically by Oster et al. [12j to model 
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waveguide arrays in a nonlinear photonic crystal. Earlier, the same equation was ob- 
tained by Claude et al. [6] for modelling of slowly varying discrete breathers in the 
Fermi-Pasta-Ulam lattices using asymptotic multi-scale expansions. Very recently, the 
generalized DNLS equation was rederived again by Abdullaev et al. [1] in a more specific 
context of the GP (Gross-Pitaevskii) equation with a periodic potential and a periodic 
sign- varying nonlinearity coefficient. If the coefficient in front of the onsite cubic non- 
linear term of the DNLS equation vanishes thanks to the sign-varying nonlinearity (that 
is, if /3 = in (fT|)), the authors of [1] incorporated other intersite cubic nonlinear terms 
from a decomposition involving Wannier functions. 

In what follows, we focus on the specific applications of the generalized DNLS equa- 
tion ([I]) in the context of Bose-Einstein condensates in a lattice. Therefore, we consider 
the GP equation with a periodic linear potential and a periodic sign-varying nonlinearity 
coefficient in the form 

id t y = -d^ + V{x)V + G{x)\y\ 2 ^, (2) 

where V(x) and G{x) are smooth, 27r-periodic functions on M. To make all arguments 
precise, we assume that 

V(-x) = V(x), G(-x) = -G{x), xel. (3) 

In this case, (3 = and our main result states that the intersite cubic nonlinear terms in 
the generalized DNLS equation ([1]) appear beyond the applicability of the DNLS equa- 
tion in the tight-binding approximation and hence must be dropped from the leading 
order of the asymptotic equation. Instead of these terms, the onsite quintic nonlinear 
term must be taken into account to balance the linear dispersion term in the quintic 
DNLS equation 

itn = a(c n+ i + Cn-i) + x\cn\ 4 c n , (4) 

where (a, x) are constant parameters which can be computed from analysis of the GP 
equation ((2|) with potentials (J3J) - 

Note that the approach leading to the DNLS equation is general and can be applied 
to other 27r-periodic functions V(x) and G(x). In a general case, (3 / and the onsite 
cubic nonlinear term is the only nonlinear term, which must be accounted in the cubic 
DNLS equation 

ic n = a(c n+ i + c n _i) + (3\c n \ 2 c n , (5) 

at the leading order of the asymptotic expansions. 

To compare the outcomes of the generalized DNLS equation (JTJ) with those of the 
quintic DNLS equation Q, we study bifurcations of gap solitons in the semi-infinite 
band gap. We show analytically that a and \ has equal negative signs in the semi- 
infinite band gap so that the quintic DNLS equation ([3]) always has a ground state, 
indicating that bifurcation of a gap soliton in the semi-infinite gap always occurs in the 
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GP equation ([2]) with potentials (J3j) - Recall that this bifurcation does not occur if the 
nonlinearity coefficient is sign-definite and positive, see Pankov |13j . 

In contrary to the predictions of the quintic DNLS equation ^ , we also show that 
the corresponding version of the generalized DNLS equation ([I]) does not admit localized 
solutions for any values of a and 7 (when [3 = 5 = 0) at least in the slowly varying 
approximation. A numerical test with particular potentials 

V(x) = Vq{1 - cos(x)), G{x) = G sin(x), (6) 

indicates that the gap solitons do exist in the semi-infinite gap independently of the 
signs of Vo and Go- The rigorous proof of existence of localized solutions in the GP 
equation ([2]) with potentials ([3j) in the semi-infinite band gap is beyond the scopes of 
this work and is a subject of an ongoing work [15J . 

We also inspect another asymptotic reduction of the GP equation (J2]) to the con- 
tinuous nonlinear Schrodinger (CNLS) equation, see review of asymptotic reductions of 
the GP equation with a periodic potential in Pelinovsky [16] . We show that the cor- 
responding CNLS equation also has a focusing quintic nonlinear term, which supports 
the same conclusion on bifurcation of a gap soliton in the semi-infinite gap. 

We note that reductions to the DNLS and CNLS equations were recently justified 
with rigorous analysis both in the stationary and time-dependent cases, see works [17]ll8j 
in the context of the DNLS equation and works [U [7] in the context of the CNLS 
equation. Therefore, it is a matter of a routine technique to formalize arguments of our 
paper. 

We shall add that the literature on the GP equation ([2]) is rapidly growing in physics 
literature. The GP equation with a periodic nonlinearity coefficient was considered by 
Fibich at al. [8], where no linear potential V(x) was included and the mean value of 
G(x) was assumed to be nonzero. For the same equation, Sakaguchi & Malomed [20] 
derived a quintic CNLS equation in a slowly varying approximation of a broad soliton. 

A more general equation with both linear and nonlinear periodic coefficients was 
studied by Bludov et al. in [3J H] , where gap solitons were approximated numerically. 
It was shown in these works that bifurcations of small-amplitude gap solitons near the 
lowest band edge depend on the sign of the cubic coefficient in the effective CNLS 
equation. Using perturbation theory, Rapti et al. [H] studied existence and stability of 
gap solitons in the semi-infinite gap for the GP equation with small linear and nonlinear 
periodic coefficients. 

The paper is organized as follows. Section 2 justifies the asymptotic reduction of the 
Gross-Pitaevskii equation ([2]) to the quintic DNLS equation (jlj). Section 3 gives results 
on existence of stationary localized modes in the generalized DNLS equations (pQ) and 
([4]) and discusses the relevance of previous works [22] and [I]. Section 4 justifies the 
asymptotic reduction to the quintic CNLS equation. 

Acknowledgements: J. B.-B. has been partially supported by grants PCI08-0093 
(Consejeria de Educacion y Ciencia de la Junta de Comunidades de Castilla-La Mancha, 
Spain), PRINCET and FIS2006-04190 (Ministerio de Educacion y Ciencia, Spain). J. 
B.-B. also would like to thank the Mathematics Department at McMaster University 
for their hospitality during his visit there. 
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2 Reductions to the DNLS equation 

To consider the tight-binding approximation and reductions to the DNLS equation, we 
assume that 

V(x)=e~ 2 Vo(x), (7) 

where e is a small parameter and Vq is a smooth, 2-7r-periodic, and even function on R. 
In what follows and without loss of generality, we set 

Vb(0) = and Vo'(0) = 2, (8) 

so that Vq(x) = x 2 + C(x 4 ) as x — ► 0. For particular explicit computations, we consider 
the standard example 

V Q (x) = 2(1 - cos(x)) = 4 sin 2 (|) . 

The limit e — > is generally referred to as the semi-classical limit |10j . 
Let k) be the Bloch function of 

L^(x;k) = E(k)^(x;k), L = —d 2 + e~ 2 Vo(x), (9) 

for the lowest energy band function E(k). It is known (see review in Pelinovsky et al. 
[TTj) that E(k) and ^(x,k) satisfy 

E(k) = E{k + 1) = E(—k), k£R 

and 

<S(x;k) = e- 2nki ^(x + 2ir;k) = $>(x;k + 1) = ®{x;-k), x£R, k G R, 
so that one can define the Fourier series decompositions 

E(k) = E n e 2 ™ k \ *(x; k) = £ ^(x)e 2 ™ fc \ 

with real-valued Fourier coefficients satisfying the reduction 

E n = E_ n , ip n (x) = 4> (x - 27rn), n G Z. 

Functions {f/Vt(£)}nez are referred to as the Wannier functions. For the lowest energy 
band, these functions form an orthonormal basis in a subspace of L 2 (R) associated with 
the lowest energy band, enjoy an exponential decay to zero as |x| — > oo and satisfy the 
system of differential equations 

(l - £ ) Mx) =J2^n (Mx) + 4>-n{x)) , xeR. (10) 

n>l 
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Thanks to orthogonality and normalization of the Wannier functions, we infer that E n 
can be computed from the overlapping integrals 



E n = (L^o^n 



$b( x )i>'n( x ) + e 2v o(x)ipo(x)tp n (x) 



dx, n E N. 



:n) 



For the semi-infinite gap and for even potentials, Wannier functions {ip n (x)} n £i are 
strictly positive and even on R. It is proved with the standard technique in the semi- 
classical limit e — > (see review in Aftalion & Helffer [2]) that the Wannier function 
ipo{x) can be approximated near x = by the normalized Gaussian eigenfunction of 

„2- 



-c^ + ^-j M x ) = -tpo(x), ie 



or explicitly, 



M x ) 



[ire 



)i/4 



e 2e x E 



This approximation suggests that 

En ~ — 



tpo(x) ~ V'o(^)) near x = 0, 



(12) 



(13) 



where we have used the notation A{e) ~ -B(e) for two functions of e near e = to 
indicate that A{e)/B{e) — > 1 as e — > 0. To obtain approximations for the overlapping 
integrals (jlip . one need to proceed with the WKB solution 



Mx) ~ A(x)e-^ S(a '' )cte ', x E (0,2vr), 



(14) 



where 



1 

(vre) 1 ^ 



cxp 



, xE(0,2tt). 



/o 25(a/) 

The WKB solution (|14p is derived by neglecting the term A" (a;) in the left-hand-side of 
(|10P and by dropping the right-hand-side of (|10|) thanks to the hierarchy of overlapping 
integrals in 

... < \E 2 \ < |ii| < |£ |- (15) 

In addition, to derive the explicit expression for A(x) we have replaced Eq by 1/e and 
used the matching condition of ipo(x) with il>o(x) as x j 0. Note that the expression for 
A{x) diverges as x f 27r. 

Thanks to the explicit formulas and the symmetry of ipo(x) on R, the first overlapping 
integral is computed as follows 



E-i 



Mx) -dl + e~ 2 Vo(x) - E ) ^ (x - ^)dx 



4^(7r)^(7r) + 2 / Mx - 2vr) -Sg + e^Or) - S U (x)dx. 



5 



Neglecting the integral (thanks again to smallness of the right-hand-side of (|10p on 
(— oo, — 7r]) and substituting the WKB solution ([H]) at x = it, we infer that the leading 
order of the first overlapping integral is given by 

E x ~ #o(7r)^(7r) = exp (- * V^)dz + jT ■ ( 16 ) 

For instance, if Vq{x) = 2(1 — cos (a;)), then 

5(x) = 2sin(f), A(x)= 1 xG(0,27t), 

v2/ (vre) 1 / 4 cos (|J 

so that 

A 16 -S 

ai 7rl/2 e 3/2 e • 

Similarly, one can establish the hierarchy of other overlapping integrals in (|15|) . See 
Helffer [10] for rigorous justification of the WKB solutions above. 

To deal with nonlinear terms, we compute the integral involving G(x)iPq(x) as e — » 0. 
This integral can be computed with the use of the Gaussian approximation (|12|) - (|13j) . 
thanks to the fast decay of ipo(x) on R and the smoothness of G(x) on R: 

/ G(x)^(x)dx ~ 1 / G(x) e -^dx ~ —A— G(0). (17) 
Tre Jr (27re) i /- ! 

The overlapping integrals involving homogeneous quartic powers of ipo{x), ipo(x — 2ir), 
etc. are much smaller compared to the integral (|17j) . thanks again to the fast decay of 
-00 0*0 on R. 

2.1 Reduction to the cubic DNLS equation 

Asymptotic reduction to the cubic DNLS equation holds for G(0) ^ 0. Computations 
(|16p and ([TT]) suggest the use of the scaling transformation 



with a new small parameter 

^ = _l_ e -f/;^M^, (is) 

for asymptotic solutions of the Gross-Pitaevskii equation ([2]). To give main details, let 
T = fit be slow time and decompose 



^o = ^c n (T)^„(x), 
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for some coefficients {c n } ng z to be defined. The remainder term \&i satisfies 



id&i = (L — Eo)Vi + ( ~ id ™ + A 4 1 S ^ m ( Cn+m 

nez \ men 



+ Cr, 



V /2 G(x)|*o + M*i| 2 (*o + /i*i)- 

Coefficients {c n } n< =z are uniquely denned by the orthogonality condition 

{4) n ,^i) = for all n G Z, 

which ensures that is in the orthogonal complement of the subspace of L 2 (M) corre- 
sponding to the lowest spectral band of operator L. Orthogonal projections to {VVijnez 
truncated at the leading-order terms as jx — > take the form of the cubic DNLS equation 



where 



a 



a(c n+1 + Cn-x) + (3\c n \ 2 c n , 

™ 1 - S'(x) 



(j, l E\ ~ — A^/Vq{tt) exp 



S(x) 



(27T) 1 / 



T7^G(0), 



(19) 

(20) 
(21) 



thanks to the fact that other overlapping integrals in the linear and cubic terms are 
smaller. Rigorous justification of the cubic DNLS equation (|19|) on a finite time interval 
is proved by Pelinovsky & Schneider [18J, where the main result is formulated in space 
H X (K) defined by the norm ||#|| w i( R ) := y/{(L + J)*, where L = -<9 2 + F(x). 

Theorem 1 Assume that V{x) is given by TJJ), G(0) ^ 0, and \x is given by U8\). Let 
{c n (T)} n6 2 G C 1 (M, / X (Z)) 6e a global solution of the cubic DNLS equation Iil9\) with 
initial data {c n (0)} n6 z G Zp(Z) /or any p > 5. Let ^0 S ^(R) satisfy the bound 



* _ e i/y/2^ Cn(0) ^ r 



1/4. ,3/2 



n 1 



for some Co > 0. There exists hq > 0, To > 0, and C > 0, such that for any 
fi G (0,/io)j i/ie CP equation |1|) u>ii/z initial data ^(O) = ^0 a solution ^(t) G 
C([0, To//i],W 1 (]R)) satisfying the bound 



VtG [0,T //i] : 



VraGZ / 



< Ce 1 / 4 /, 372 - 



h 1 
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2.2 Reduction to the quintic DNLS equation 

If G(0) = 0, and we assume that G(x) is odd on R, then (3 = and the DNLS equation 
(|19p becomes a linear equation. A modified asymptotic solution is needed to incorporate 
the leading order of the asymptotic expansion. We will show that the modified scaling 

tf(s,t) = e^V 74 f* + Ai 1/2 ^i + ^2) e^ ot , 



will reduce the GP equation ((2J) to the quintic DNLS equation (jij) if G'(0) / 0. Again, 
let T = lithe the slow time and define \&o arm ^1 by 

^0 = Cn{T)Mx), *1 = ^ \Cn(.T)\ 2 C n (T)<p n (x), 

where (p n {x) = 0o(x — 27rn), n £ Z is a solution of 

(L - Eq)<Po(x) = -e-V 2 G(x)i%(x), xeR, (22) 
under the orthogonality condition 

f G(x)$(x)dx = 0. (23) 

The remainder term ^2 satisfies 



(L - £ )*2 + ( ~ id n + V 1 Yj 



i<9t*2 = (L- E )^2 + / A -iCn + A* S m (c„+ m + C n _ m ) 



E^(i c ™i 2c «)^+^ v v i/2 G(x) 



X Mtfo + M 1/2 ^l + ^2| 2 (^0 + ^ 1/2 ^l + M*2) - Y M^^n J • 

Orthogonal projections to {ipn}nez truncated at the leading-order terms as fi — > result 
in the quintic DNLS equation ([J]) with the same expression for a as in (|20p and the 
following expression for 

X = 3e- 1 / 2 f G(x)i,l{x)w{x)dx. (24) 

The justification of the quintic DNLS equation dH) relies on the two facts. 

Lemma 1 Under condition there exists a solution <po(x) of the inhomogeneous 
equation A22\) so that \ is bounded and nonzero as e —* 0. 
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Proof. First we note that if G(x) is odd and ipo{x) is even on R, then <po(x) is odd on 
R, so that the integral in (|24p is generally non-zero. Moreover, using the inhomogeneous 
equation (|22p . we infer that 

X = -3((L- E o )0o,0o), (25) 

so that x < for the lowest band of L, since (L — Eq) is positive definite if Eq is at the 
bottom of the spectrum of L in the limit e — ► 0. To show that x is bounded as e — > 0, we 
can use again the Gaussian approximation (fT3j) and find solutions of the inhomogeneous 
equation (|22|) near x = in the form 

eV2G'(0) . . 

(p (x) r — — -Tjr-xe 2, , near x = 0. (26) 

As a result, 

x ~~w^> (27) 

and we see that x is bounded and negative as e — > 0. ■ 

Lemma 2 Under condition {2j), £/te largest overlapping integrals from the cubic term 
|^o| 2 ^o, 

G(x)^g(a?)^o(s-27r)dx, / G(x)^g(x)^(x - 2vr)dx, (28) 
are smaller than e 1 / 2 /^. 

Proof. First, we note that if G(x) is a smooth, 27r-periodic, and odd function, then 
G(x) is also odd with respect to the point x = it, so that 



G(x)%(x)%{x - 2ir)dx = 0. (29) 
To consider the other nonzero integral in (|28p . we write 

G(a;)^g(x)^o(z-27r)dx= I / +/ ) G(x)$(x)ip (x - 2w)dx. (30) 



The second integral on [it, oo) is much smaller than the first integral on (— oo, tt] thanks 
to the faster decay of i/)q{x) compared to tpo( x ) on As a result, we deal only with 
the first integral, which we rewrite as follows: 

G(x)$(x)'$o(x - 2n)dx (31) 

12 / Mx-27r)^-dl + e- 2 V o (x)-E o ^0o(x)dx 



e l/2 



M*)<Po(*) + ^(7r)^(7r)] - e 1 / 2 / £o(x) (-d 2 + e' 2 ^^) - E ) M% ~ 2vr)dx, 
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where we recall again that i/jq(x) is even on R. In view of equation (jlOp . we have 

o (x) (-<9 2 + e~ 2 V (x) - £ ) M x ~ 27T )dx 
^2 E n 0o(x) (tp n+1 (x) + -4>- n+ i(x)J dx 

J — oo ^ ' 



n>l 



~ Ei / 0o (x)tpo(x)dx = -Ei / (f (x)ipo(x)dx, 

J — OO J 7T 

where the last equality is due to the fact that i/jo(x) is even and 0o(x) is odd on M. 

Thanks to the fast decay of ipo(x) and <po(x) on R, the second term in pip becomes 
smaller than e 1 / 2 ^! = e 1 / 2 ^a, where a is given by (|20p . 

Boundary values of ipo(x), 0o( x ) an d their derivatives at x = ir in the first term in 
(|3T|) can again be computed from the WKB solutions for tpo( x ) and <^o(;r)- 

For solutions of the inhomogeneous equation ([22]) . we substitute 

<^o(x) ~ B{x)e~^o s(x ' )dx \ x e (0,27r), 

where S'(x) = v/Tfo(x) and satisfies the first-order differential equation 

e(6S(x)B'(x) + 3S'(x)£(x) - - 8S 2 (x)B(x) = -e 3/2 G(x)A 3 (x), 

where the term B"{x) is neglected and Eq ~ 1/e is used. Solving the differential equation 
with the integration factor, we obtain 

b ^ - M) exp (s r s(i,)<fa ' + s r ■ « <°' 2i) - 

with 

c{x) =-—J svw) exp L s{x )dx " 6 L -^r dx ) dx ■ 

Using the Laplace method for computing integrals, we obtain a correct behavior of <po(x) 
near x = that matches the previous calculation ([261) : 

- - 6(7re)3/4xl/ 3 e 2e y y /3e 3> )d y 

e 1/2 G'(0) 

~ ; TTTT X ^ 2e ) near X = V. 

8(^e) 3 / 4 

As a result, we have 

e 1/2 ^(tO^oO) +^o( 7r ) ( ^o(7r) 
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where 



G(x)A 3 (x) f 4 [* , 1 f t- 5V) M , 



For instance, if V(x) = 2(1 — cos(aj)), we obtain 



e 1/2 VA)0)$)(» +^o( 7r ) ( £o(7r) 
which is clearly smaller than 



4 _48 



G(x) 



9vre y sin 2 / 3 (f ) cos 10 / 3 (f ) 



e-£ cos ^dx, 



7rV4 e 



This completes the proof of Lemma [2j ■ 

Using the approach from Pelinovsky Sz Schneider in [18], we can justify the quintic 
DNLS equation (JU) on a finite time interval, according to the following statement. 

Theorem 2 Assume that V(x) and G{x) are given by and ([?[), G"(0) 7^ 0, and is 
given by i tiffj) . Ze£ {c n (T)} ng ^ G C 1 (M, Z (Z)) 6e a global solution of the quintic DNLS 
equation initial data {c n (0)} ne ^ £ Z 2 (Z) /or any p > ^. Let \fo £ satisfy 

the bound 



w 1 



/or some Co > 0. There exists ^0 > 0, T$ > 0, and C > ; suc/i t/iai for any 
[i S (0, [io), the GP equation (0j with initial data ^(0) = has a solution ^(t) G 
C([0, To/^],W 1 (M)) satisfying the bound 



Vie [0, Tq/ax] : 



< Ce-^V 374 - 



H 1 



Remark 1 Note that the results of Theorems [7] and also hold for the piecewise- 
constant Kronig-Pennig potential Vq{x) after minor modifications required because of 
a different algebraic factor of e in the definition of \jl 

3 Localized solutions of reduced equations 

We shall consider the stationary solutions of the quintic DNLS equation (J3|), where 
the coefficients a and \ are computed asymptotically for the lowest energy band of 
L = -d 2 x + e- 2 Vb(x) in the semi-classical limit e —* 0. 
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Let Cn(T) = 4> n e tQT for a real parameter and a real-valued sequence {4>n}n& 
and obtain the stationary quintic DNLS equation 

a(> n+ i + n _i) + = fi&i, n G Z. (32) 

The hierarchy of overlapping integrals implies that the energy band function E(k) 
is given at the leading order by 

E(k) ~ E Q + 2^i cos(2vrA;) + ... 

Since = is the minimal point of E(k) for the lowest energy band, we have E\ < 
so that a < 0. See also ([20]) . where q < is computed in the limit e — » 0. On the other 
hand, representation (|25p implies that x < f° r the lowest band. See also (|27p . where 
X < is computed for Cr'(0) / as e -> 0. The semi-infinite gap corresponds to the 
interval Vt <2a. 

Localized solutions of the stationary quintic DNLS equation (|32[) can be obtained 
from a minimization of the energy functional 

ngZ 

subject to a fixed N = X^nez^n- According to Theorem 2.1 of Weinstein [23], there 
exists a minimizer of H (called a ground state) for sign(a) = sign(x) with sign(S7— 2a) = 
sign (a). Monotonic exponential decay of the sequence {^ n }nez to zero as in ±oo was 
shown in Theorem 1.1 of Pankov p3] (where the cubic DNLS equation was considered 
without loss of generality). Note that the localized solution also exists if sign(a) = 
— sign(x) for sign(f2 + 2a) = — sign(a) thanks to the staggering transformation 

<\>n (-l) n <pn, X^-X, 0->--Q, a -> a, 

that leaves solutions of (|32p invariant. Therefore, the localized solution is not monoton- 
ically decaying if sign(a) = — sign(x). For the semi-infinite gap, we have shown above 
that a and x have equal negative sign, so that a localized solution of the stationary 
quintic DNLS equation f|32|) exists in the semi-infinite gap for £1 < 2a. 
Consider the stationary GP equation 

-&'(x) + V(x)<5>(x) + G(x)$ 3 (x) =£j$(x), x£R, (33) 

which is derived from the GP equation ([2]) from fy(x,t) = <J?(x)e _ia;< . Persistence anal- 
ysis of gap solitons in Pelinovsky et al. [17] gives the following result. 

Theorem 3 Let V{x) and G(x) satisfy {3|) and ^ and G"(0) ^ 0, and [i is given by 
fl8\) . Let {4>n}nei G I (Z) be a ground state of the stationary quintic DNLS equation 
nnfy for Q < 2a. There exists no > and C > 0, such that for any fi G (0,//o)> the 
stationary GP equation \33\) with lu = Eq + fj,Q has a solution <3? G Ti. 1 (M.) satisfying the 
bound 

< Ce" 1 V /4 - 

Moreover, <j)(x) decays to zero exponentially fast as \x\ — > oo. 
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Remark 2 One can also prove existence of gap solitons in the semi-infinite gap of the 
GP equation Hjj with potentials (Qp in the opposite limit of large- amplitude gap solitons 
using the Lyapunov- Schmidt reduction method. See Sivan et al. ]21^ for an example of 
this technique for the GP equation with a periodic linear potential V(x) and a constant 
nonlinearity coefficient. 

To summarize, from Theorem [3j we predict existence of gap solitons in the semi- 
infinite gap for any even V(x) and odd G(x) with G"(0) 7^ 0. To illustrate the existence 
numerically, we solve the GP equation by using the so-called imaginary time method 
[9]. As approximations of localized solutions evolve along the imaginary time, iterations 
converge to the ground state of the stationary GP equation ([33]) . 

We have developed a Fourier pseudospectral scheme for the discretization of the 
spatial derivatives combined with a split-step scheme for iterations in the imaginary 
time, see implementation of this method by Montesinos & Perez-Garcia [llj. In other 
words, solutions of 

d t U{x,t) = (A + B)U, 

with 

A = -d xx , B = V(x) + G(x)\U\ 2 , 

are approximated from exact solutions of the problems dtU = AU and dtU = BU. By 
using the symmetric (second-order) split-step method, whose equation is 

U{x, t + r) = e TA/2 e TB e rA/2 U{x, t) + 0(r 3 ), (34) 

we calculate a localized solution of the stationary GP equation ([33]) as t — » 00. Figure 
[T] shows the branch of gap solitons bifurcating to the semi-infinite gap (left) and a 
particular profile of the localized solution (right) that corresponds to the point on the 
solution branch on the left. 

We note that the numerical scheme we have used has many advantages. First, it is 
more accurate than finite-difference numerical methods. Second, the Fourier transform 
can be computed by using the fast Fourier transform. Finally, the L 2 -norm of the 
localized solutions is preserved during the time iterations so that the L 2 -norm of a gap 
soliton along the solution branch can be fixed by the starting approximation. 

In the end, we note that existence of localized solutions in the stationary GP equation 
([33]) for any smooth 2-7r-periodic even V(x) and odd G(x) in the semi-infinite gap of L 
can be proved using the variational theory by a modification of arguments in [13]. This 
modification is a subject of an ongoing work [15]. Numerical evidences of existence 
of gap solitons in the semi-infinite gap for sign-varying nonlinearity coefficients can be 
found in [3l|4]. 

3.1 Comparison with the generalized DNLS equations 

Let us compare our main conclusion with the prediction of the stationary generalized 
DNLS equation considered by Abdullaev et al. \X\. For the case of odd nonlinearity 
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Figure 1: The solution family of gap solitons for Go = —10 and Vq = 6 in ©: The 
L 2 -norm N versus u> (left) and the spatial profile of gap soliton corresponding to marked 
point with a black circle (right). 

coefficient G(x), this stationary equation is written in the form 

Oi((j) n+ i + <f> n -l) +7(3<^(0 n+ i - <t> n -l) ~ 4>n+X + 4>n-l) = ^f>n, n £ Z, (35) 

where a is the same as in (132|) and 7 is proportional to the overlapping integral (1301) . 
Note that the cubic term in (|35p is slightly different from the one in ([T|), which holds 
for even nonlinearity coefficient G{x) pQ. We also note that [3 = 5 = thanks to (123ft 
and (I29p . The energy functional is now written as follows: 

H = ^ (<X<t> n <l>n+\ + l^ni^n+l ~ <t>n-\)) ■ 

While we are not able to prove that the stationary DNLS equation (|35p admits no 
localized solutions for any signs of a and 7, we can simplify the problem in the slowly 
varying approximation, which is also referred to as the continuum limit of the lattice 
equation. To this end, we assume that the following expansion makes sense 

0n±l = <f>(x n ) ± H'{x n ) + hi 2 4>"{x n ) + 0(h 3 ), 

where x n = hn, n £ Z and apply the scaling 

a = 2ha, tt - 2a = 2h 3 tt. 

At the leading order, the difference equation (|35j) becomes the second-order differential 
equation 

a<j)"(x) - 70' (s) [(<^'(x)) 2 + 2>4>{x)<t>" {x)] = fu/>(x), xeR, (36) 
which has the first integral 

/ = ^a(0'(x)) 2 - 7 <A(x)(0'(x)) 3 - l -£l<p\x). 
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We note that 1 = for localized solutions and that no turning point Xo £ R with 
4>(xq) > and 4>'{xq) = exists. As a result, the trajectory departing from the critical 
point (4>,4>') = (0,0) in the first quadrant of (0, </>') remains in the first quadrant and 
goes to infinity. As a result, no classical localized solutions of the differential equation 
(|3f)j) exist. Thus, we have the following result. 

Proposition 1 Stationary generalized DNLS equation \35\) for any coefficients a, 7, 
and Q admits no localized solutions in the slowly varying approximation. 

We conclude that the stationary generalized DNLS equation ([TJ gives the opposite 
(wrong) conclusion to the bifurcation problem of localized solutions in the semi-infinite 
gap, compared to the stationary quintic DNLS equation (j32l . 

It is even more problematic how to interpret the modification of the generalized 
DNLS equation (pQ) by Smerzi & Trombettoni [22], where the onsite cubic nonlinear 
term P\c n \ 2 c n was replaced by (5\c n \ 2p c n with p < 2. If f3 ^ 0, the justification of the 
cubic DNLS equation (119D in Theorem Q] leaves no hope to have p < 2 in the generalized 
DNLS equation and to account the intersite cubic nonlinear terms at the same order 
as the onsite cubic nonlinear terms. Thus, we have to conclude that the generalized 
DNLS equations considered in [1] and [22] (and implicitly in [6] and [12]) are invalid for 
potential V(x) in ([7]) in the tight-binding approximation as e —* 0. 

4 Reductions to the CNLS equation 

Let us now consider the potential V(x) in the GP equation ([2]) without assumption 0. 
Spectral bands are generally of a finite size, so that we can simplify the GP equation ([2]) 
if the bound state has small amplitude near the band edge. This asymptotic reduction 
leads to the continuous nonlinear Schrodinger (CNLS) equation justified by Busch et 



To give main details, let Eq be the lowest band edge of operator L = —d 2 . + V(x) 
corresponding to the 27r-periodic L 2 -normalized eigenfunction ^0 £ ^per(0>27r). Since 
the second solution of L^> = Eq^ is linearly growing, the subspace Ker(L — EqI) C 
Lp er (0,27r) is one-dimensional. Looking at the Fredholm alternative condition for the 
inhomogeneous equation 



we infer that there exists a unique 27r-periodic function ^1 £ Lp er (0, 2ir) in the orthog- 
onal complement of Ker(L — EqI). If V(x) is even on R, then ^0(2^) is even and ^i(x) 
is odd on R. In addition, if G(x) is an odd 27r-periodic function, there exists a unique 
odd 2-7r-periodic solution of the inhomogeneous equation 



al. 0. 



- &[(x) + V(x)^i(x) - E ®i(x) = 2%{x) 



(37) 



V 2 '(x) + V(x)V 2 (x) - E V 2 (x) 



G(x)^ (x) 



(38) 
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that also lies in the orthogonal complement of Ker {L — EqI) . Equipped with these facts, 
we are looking for an asymptotic solution of the GP equation ([2]) using the decomposition 

*(x,t) = e 1/2 {A(X,T)^ (x)+e{A x (X,T)^ 1 (x) + \A(X,T)\ 2 A(X,T)^ 2 (x)) 
+s 2 ^(x,t)) e~ iEot , 



where e is a small parameter, X = ex and T = e 2 t are slow variables, and ^(x,t) 
satisfies the time evolution equation 

idt^f = (L - E )^ - iA T ^o - e (Axt^i + {\A\ 2 A) T ^ 2 ) 

-A xx (* + 2^i) - 2(\A\ 2 A) X % - e (Axxx^i + (\A\ 2 A) XX * 2 ) 

+G(x)e- 1 (\AV + e (A x Vi + \A\ 2 A* 2 ) + e 2 ^\ 2 
(At> + e (A X V 1 + \A\ 2 A* 2 ) + e 2 §) - |,4| 2 ^' 



Projecting the right-hand side to an d truncating at the leading-order terms, we 
obtain the CNLS equation 



iA T = aA xx + x \A\ A A + 1 (\A\ 2 A) x , (39) 



where 



/•27T 

a = -1-2 / ^/[(x)^ (x)dx, 
Jo 



2n 



X = 3 / G(x)V 6 (x)t> 2 (x)dx, 
Jo 

7 = -2 1 ty' 2 (x)$>o(x)dx + / G(x)^l{x)^x(x)dx. 
Jo " Jo 

Justification of the generalized CNLS equation (|39p can be developed similarly to the 
work of Busch et al. [5 J . While it may seem that the generalized CNLS equation (|39p 
contains both the quintic and the cubic derivative terms, we obtain that 

r-2-K 

7 = / (-2^ 2 (x)^ (x) + G(x)^l(x)^ 1 (x))dx 
Jo 

/■27T 

= -/ (2y' 2 (x)y (x) + y 1 (x)(-d 2 x + V(x)-E )y 2 ( y x))dx 
Jo 

= - / (2* 2 (x)V' (x) + V 2 (x) {-d 2 x + V(x) - E ) *i(x)) dx = 0. 
Jo 

Therefore, the generalized CNLS equation (|39p is just the quintic CNLS equation 

iA T = aA xx + x \A\ 4 A. (40) 
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For stationary solutions with A(X,T) = a(X)e where Q, and a(X) are real- 
valued, we obtain the stationary quintic NLS equation in the form 

aa"(X) + xa 5 ( x ) = Va(X), IeK. (41) 

Similarly to the case in the tight-binding approximation, we note that a < and % < 
for the semi-infinite gap since a = —^E"(0) < 0, where E(k) is the energy band function 
for the lowest energy band, and 

X = 3 / G(x)^l(x)ty 2 (x)dx = -3 / tf 3 (a;) {-d 2 x + V(x) - E Q ) ^ 2 {x)dx < 0. 
Jo Jo 

The stationary quintic NLS equation (|4ip has a positive definite soliton for sign(a) = 
sign(x) with sign(S7) = sign(a), that is for O < in the semi-infinite gap. 
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